Phase-separation transition in liquid mixtures near curved charged objects 
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We study the thermodynamic behavior of nonpolar liquid mixtures in the vicinity of curved 
charged objects, such as electrodes or charged colloids. There is a critical value of charge (or 
potential), above which a phase-separation transition occurs, and the interface between high- and 
low-dielectric constant components becomes sharp. Analytical and numerical composition profile 
are given, and the equilibrium front location as a function of charge or voltage is found. We further 
employ a simple Cahn-Hilliard type equation to study the dynamics of phase-separation in spatially 
nonuniform electric fields. We find an exponential temporal relaxation of the demixing front location. 
We give the dependence of the steady-state location and characteristic time on the charge, mixture 
composition and ambient temperature. 
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Situations where charged objects, such as electrodes 
or colloids, are found in liquid environments are abun- 
dant in science and technology. In ionic mixtures, ex- 
periments and theory show that ions migrate toward 
the object and lead to screening of the electric field. In 
nonpolar liquids, the situation is different: the decay of 
electric field far from the objects depends on the geom- 
etry of all conducting surfaces, and may be long range. 
When a nonpolar liquid mixture is under the influence 
of a uniform electric field E, the theories of Landau [l[ 
and later Onuki [2] showed that the critical temperature 
can change by a small amount, proportional to E 2 . Ex- 
periment by Debye and Kleboth [3] partially confirmed 
the theory. 

However, here we show that the situation in spatially 
nonuniform electric fields, occurring when liquid mix- 
tures are found under the influence of curved charged 
surfaces, is quite different. When the temperature T 
is larger than the critical temperature T c , the mix- 
ture exhibits smooth composition variations. This di- 
electrophoretic behavior is reminiscence of the effect of 
gravity. Q For a homogeneous mixture below T c , there 
are two scenarios: if the charge density is small, there 
are still weak composition gradients. On the other 
hand, large enough charge leads to a phase-separation 
transition, where the liquid with high-dielectric con- 
stant is close to the high field region while the liquid 
with low dielectric constant is pulled away, and the co- 
existing domains are separated by a sharp composition 
front. 

The phase transition described below occurs in 
systems described by bistable free-energy functionals 
giving rise to a phase diagram in the composition- 
temperature plane divided into two regions: homoge- 
neous mixture and a phase-separated state. In order 
to be specific and to facilitate the connection to exper- 
iment, we consider the following binary mixture free- 
energy density f m = kTf m /Nv , where 



f m = 01og(</>) + (l 



dog(l 



+ N X cf>(l-cb) (1) 



This symmetric (Na — Nb = N) free energy is given in 
terms of the A-component composition (p (0 < (j> < 1) 
in a mixture of A/B liquids, and the so-called Flory 



parameter x ~ l/^-H Simple liquids have N = 1, 
while polymers have N > 1 monomers, each of vol- 
ume vq. k is the Boltzmann constant. The critical 
point is given at (cj) c ,(Nx)c) = (1/2,2). In the ab- 
sence of electric field, the mixture is homogeneous if 
T > Tt, and unstable otherwise. The transition (bin- 
odal) temperature T t at a given composition is given 
by T t = (A^x)cT c [log(0/(l-0))/(20-l)]- 1 .[5] The 
phase transition does not depend on the exact form of 
f m , and appears in a Landau series expansion of Eq. 
([T]) around <p c , or in any other similar "double- well" 
free-energy functional. 

As is shown below, the effect of electric fields is large 
if they originate from curved charged surfaces. In this 
work we consider for simplicity surfaces with fixed cur- 
vature: a charged spherical colloid, a charged wire or 
two concentric cylinders, and the "wedge" capacitor, 
made up from two flat and nonparallel surfaces. Fixed 
charges on the conductors, fixed potentials, or a combi- 
nation of the two are considered by us. When the mix- 
ture is in the vicinity of a charged object with a fixed 
surface charge, the total dimensionless free-energy is 
/ = fm+U, where / os = {Nv /kT)[{l/2)e^){V^f\ is 
the dimensionless electrostatic free energy density. [l|, Q 
Note that we do not include any direct short- or long- 
range interactions between the liquid and the confining 
walls. 

The equilibrium state is a solution of the two cou- 
pled nonlinear equations: 8f/5<p — and 5f/8ip = 0, 
where ij) is the electrostatic potential obeying the proper 
boundary conditions. 0, 0] The equation Sf/S'tjj = 
leads to Laplace's equation: V • (e(^)W) = 0, and 
is readily solved by the use of Gauss' law for sys- 
tems with prescribed charges on the confining conduc- 
tors and in azimuthal or spherical symmetries. For 
example, for a mixture confined between two infi- 
nite concentric cylinders of radii R\ and i?2 > Ri, 
we find E = — Vi/j = \/(2ire((f>)r), where A is the 
charge per unit length on the inner cylinder and r is 
the distance from the cylinder's center. Subsequently, 
5f„/8<f> = -~N 2 vox/(kT c )(\/(4irsr)) 2 ds/d(l). Similarly, 
E = Q I ' (A-Ke(4>)r 2 ) for a spherical colloid of radius R\ 
and charge Q, and r is the distance from the colloid's 
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center, and E = V/((3r) for a wedge consisting of two 
flat conductors with potential difference V and opening 
angle (3 between them, and r is the distance from the 
imaginary meeting point of the conductors. 

We thus arrive at a considerable simplification of the 
problem, since the expression for E obtained above al- 
lows to write a single dimensionless governing equation 
for all three cases with radial or azimuthal symmetry: 

log (^) + N X (1 - 20) - N X M^r-« - M = 

In the above, M is the dimensionless ratio between 
the maximum electrostatic energy stored in a molec- 
ular volume and the thermal energy. M is M c = 
A 2 Nvo / (16ir 2 kT c R 2 eo) for two concentric cylinders, M 
is M s = Q 2 Nv Q /(6ATr 2 kT c s a Rt) for a spherical col- 
loid, and M is e 2 M w for the wedge, where M w = 
V 2 Nv Q e Q /(4p 2 kT c Rl), V is the voltage between the 
wedge plates, and R\ is the smallest distance from 
the conductors' edge to their imaginary meeting point, 
f = r/Ri is the scaled distance from the center of the 
sphere or the inner cylinder, and e = e/eo, where £o 
is the vacuum permittivity. Finally, n is the expo- 
nent characterizing the fall of E 2 : n = 2 for concentric 
cylinders and the wedge geometries, and n = 4 for the 
sphere. The importance of curvature is exemplified by 
the appearance of i?i in the expressions for the M's. \i 
is a Lagrange multiplier needed to conserve the average 
mixture composition: (0(r)} = </>o, and 0o is the aver- 
age composition. In the case of an open system coupled 
to a particle reservoir at r — > oo, /i is the reservoir's 
chemical potential. The phase-transition described be- 
low is from a homogeneous (mixed) to a demixed state, 
and therefore it is assumed that 4>o is outsides of the 
binodal curve, namely, T > T t . 
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FIG. 1: Equilibrium profiles </>(f) for concentric cylinders 
with different dimensionless charge M c . Dashed line: T = 
0.991T C , and M c = 0.008 too small for phase separation. 
Circles: same T, but M c = 0.04. Solid line: same T, but 
M c = 0.08. Dash-dot line: M c = 0.08, but T = 0.994T C is 
a higher temperature. We took fixed average composition 
4>o — 0.4. In this and other figures, -Ri = 1, R2 = 5, e a = 
5eo, and eh = 3eq- 



Equation ^ expresses implicitly the composition 
profile <j){f). Above T c , <f>(r) has only smooth variations, 
irrespective of the value of M. Below T c [equivalently 
N\ > (Nx)c]i if M is sufficiently small, the profile <j>{r) 
is smooth, with high-0 values at small f's and low values 
at larger radii. However, there is a critical value of M, 
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FIG. 2: (a) Filled symbols: dimensionless critical charge M s * 
for phase-separation near an isolated spherical colloid as a 
function of temperature. The colloid is coupled to a reservoir 
with three compositions: (j>o = 0.2 (squares), 0.3 (circles) 
and 0.4 (crosses). Open symbols: M* for a closed cylindrical 
system with same compositions, (b) M^ vs A(j> — fit — 4>o 
from Eq. © (solid line) and from numerics (symbols). 

denoted M* , above which <p(f) exhibits a sharp jump: 
for M > M*, high- and low-0 domains coexist sepa- 
rated by a clear interface at f = R. This transition oc- 
curs generally, even when the constitutive relation e(<p) 
is linear in (f>. This is in contrast to the Landau mech- 
anism, which relies on a quadratic dependence of e on 
(f> and is responsible to a small change in T c . We there- 
fore chose the linear relation e(4>) — Sb + 0Ae, where 
Ae = e a ~Eb, and e a and Sb are the dielectric constants 
of components A and B, respectively. 

The typical demixing electric fields and surface or 
line charge density can be estimated from the values of 
M (see Figs. 1 and 2). At M s = 0.001 and using a 
molecular volume of Nv = 10~ 26 m 3 , colloid's radius 
Ri = 1 /iixi, T c ~ 300 K, e ~ 4, we find the electric field 
at the sphere's edge to be E ~ 10 6 — 10 7 V/m (surface 
potential ~ 1 — 10 V). The corresponding charge density 
is a = eE - 10~ 5 - 10~ 4 C/m 2 (total charge Q = 
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800 — 8000e). Similar values for the electric field and 
charge density appear in the concentric cylinders and 
wedge geometries. 

Figure 1 shows 0(f) for a binary mixture confined 
by two concentric cylinders for several values of the di- 
mensionless parameter M c and at two different temper- 
atures. When M c = 0.008, there is no phase separa- 
tion, and the profile is smooth. As M c increases above 
M*, phase separation occurs, and 0(f) rapidly changes 
from high to low values at the phase-separation front 
located at f = R. Further increase of M c at constant 
temperature leads to displacement of R to larger values 
and to larger composition difference between coexisting 
domains. d,Q 

Figure 2(a) shows the calculated critical value M* as 
a function of temperature for a spherical colloid coupled 
to a particle reservoir at three different compositions. 
At a given T above T t , larger values of |0 C — 0o| require 
more charge for demixing. Curves also show M* for a 
system enclosed between two concentric cylinders. No- 
tice that approaching T t , M* becomes infinitesimally 
small. For a wedge with average composition 0o close 
to the transition composition t at given temperature, 
we obtain the following approximation: 



M* 



t -0o T rf 2 / m (0 t ) 
4Ae T c dS 2 



(3) 



where x = R2/R1 and g[x) = 2(x 2 — 1 )/(x 2 — 1 — 2 lnx). 
Figure 2(b) shows M£ from this formula and compares 
it with a more accurate numerical solution. 





100 200 300 400 
time 



FIG. 3: Composition profiles <p(r , i) at several dimensionless 
times for concentric cylinders with M c = 0.32, 0o = 0.3, and 
T — 0.95T C . Inset: semilog plot of R(t). Numerical results 
(squares) and experiments of Ref. (circles, time in s). 



We now turn to describe the relaxation toward equi- 
librium. The dynamics are governed by the following 
set of equations: (TT 



n|,[i2j 



^ + u-V0 
Vu 



LV 2 5f/S(/) 






du 

~dt 



+ (u • V)u 



r/V 2 u 



(4) 

(5) 
(6) 

VP — 0V<5//<50 (7) 



u is the velocity field corresponding to hydrodynamic 
flow and 77 is the liquid viscosity. Equation ([¥]) is a conti- 
nuity equation for 0, where —LV(df /6<fi) is the diffusive 
current due to the inhomogeneities of the chemical po- 
tential, and L is the transport coefficient (assumed con- 
stant). Equation (J5J) is Laplace's equation, Eq. ^ im- 
plies incompressible flow, and Eq. Q is Navier-Stokes 
equation with a force term —6W5 f/86.\l'q 

We continue in the limit of overdamping and with the 
assumption of azimuthal symmetry. It follows that u = 
0. We use the dimensionless time t = N VoRft / '(LkT) , 
radius f = r/R\ and energy / = Nvof /kT, to ex- 
press 4> as a solution to a diffusionlike equation dcf>/dt = 
V 2 <5/ /6(f>, while satisfying Laplace's equation, where the 
"~" signs have been omitted for brevity of notation. 
The time dependence of the profile <fi(r, t), obtained nu- 
merically, is shown in Fig. 3 for several times t. 
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FIG. 4: Steady-state front location Roo vs time constant r 
for different temperatures and compositions. Stars: M c = 
0.318, 0o = 0.3 and 0.95 < T/T c < 0.99. Circles: M c = 
0.818, 0o = 0.2, and 0.89 < T/T c < 0.95. Squares: M c = 
0.818, 0o = 0.3, and 0.95 < T/T c < 0.97. Up and down 
triangles and diamond: 0o = 0.22, 0o = 0.24, and 0o = 0.26, 
respectively, M c = 0.818 and T/T c = 0.95. 

The dimensionless location of the demixing front 
changes as a function of time: R = R(t) and asymptoti- 
cally tends toward the steady-state front location R^ at 
long times. We find excellent match with an exponential 
relaxation of the form R(t) = R ao + (l — i?oo) exp(— t/r), 
as is shown in the inset of Fig. 3. 

Figure 4 shows the location of the steady-state demix- 
ing front Roc and the time constant r at several different 
values of 0o and T, and for two different values of M c . 
It is worth noting that all the points with the same M c 
seem to fall on the same line. Similarly, the dependence 
of Roo on 0o is displayed in Fig. 5(a). Clearly, the do- 
main size increases with M c at constant temperature 
and composition. Increase of 0o at constant T and M c 
increases the domain size. Figure 5(b) shows how r de- 
pends on 0o. Compositions closer to C exhibit slower 
relaxations. In addition, increase of M c leads to faster 
relaxation toward steady state. 
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It should be emphasized that this phase transition is 
not restricted to the vicinity of the critical point, and 
it occurs at all compositions, provided that the electric 
field is large enough. Moreover, field-induced prewet- 
ting could also be realized in vapor-liquid systems of 
pure substances subject to nonuniform electric fields. 
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plication of external potential, recalling that in these 
systems the electric field is inherently nonuniform. 

Lastly, we point out that the demixing transition 
creates optical interfaces, since the mixture's compo- 
nents have different refraction indices. Consequently, 
the propagation of a light beam through a mixture in a 
channel will be altered once an electric field creates op- 
tical interfaces, and this may be used to scatter, focus, 
or even guide rays in microfluidic arrays. (l9j 

We thank L. Leibler and F. Tournilhac for help in 
developing the ideas presented in this work, and D. An- 
delman for numerous useful comments. This research 
was supported by the Israel Science foundation (ISF) 
grant no. 284/05, and by the German Israeli Founda- 
tion (GIF) grant no. 2144-1636.10/2006. 



150 


o M c =0.4 

a M =0.8148 
v c 












e 100 

50 







% 


o 










0.2 



0.22 



0.24 0.26 
4>0 



0.28 



0.3 



FIG. 5: Steady-state front location Roo (a) and time con- 
stant r (b) vs <f>o for two values of dimensionless charge M c . 
Numerical solution for concentric cylinders with Ri = 1, 
R 2 = 5, andT = 0.95T c . 

There are several circumstances where the field- 
induced separation may have an important influence. 
Colloidal suspension in liquid mixtures and poly mer so- 
lutions have been extensively studied. [H, HI , [IU HI, 
17| We point out that standard wetting theory is in- 
sufficient to describe these experiments if the colloids 
are charged. The enrichment layer around the colloid is 
sensitive to the colloid's charge, and this may have an 
effect on the intercolloid interaction and hence on the 
phase behavior and the rheology of suspensions. 1(| 17 1 

A drastic change to the rheological properties is 
also predicted for a mixture confined, for example, be- 
tween two rotating coaxial cylinders (Taylor-Couette 
flow). The classic (zero field) flow profile would change 
markedly if a potential is imposed between the two 
cylinders. Once the homogeneous mixture demixes, 
most of the velocity gradient will fall on the liquid com- 
ponent with smaller viscosity. [H} A change to the lu- 
brication in microelectromechanical systems and in mi- 
crofluidic channels can be similarly brought by the ap- 
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